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ABSTRACT 

A method is presented for the analytic design of linear 
control systems. Design is accomplished by the solution of 
an algebraic set of simultaneous equations derived from the 
design specifications. Specifications considered include the 
system dam^Ang ratio, undamped natural frequency, bandwidth, 
steady state error coefficients and root sensitivity. 

The method is showr^c^^BS^fa-dl, accurate , specific and 
capable of working with any number of parameters . In 
addition, the analytic approach minimizes dependence on 
graphical techniques and trial and error solutions. 
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SYMBOLS AND ABBREVIATIONS 



CE 


Characteristic equation 


Ess 


Steady state error 


G(s) 


Forward transfer function 




Acceleration error constant 


Ka 


Accelerometer gain 


K 


Gain 


Kp or kp 


Position error constant 


Kv or 


Velocity error constant 


Kt 


Tachometer gain 


^pw 


Reasonance peak 


Superscript or 
subscript "o" 


Specified or numerical value 


w (omega) 


Angular frequency 


Wn 


Undamped natural frequency 


Wb 


Bandwidth ( -3 db) 


Wr 


Rea sonant frequency 




Polynomial in ^ 


Q 


General sensitivity parameter 


s 


Complex variable of the form cr + jw 


S 


Sensitivity 




Sensitivity of P with respect to Q 
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S^n 
q RT 

^ (zeta) 



Jr component of root sensitivity 
Wjj component of root sensitivity 

Total root sensitivity 

Damping ratio 
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1 . Introduction . 



Linear control system design can be defined as "a set of 
procedures for determining the value of system parameters 
which will cause a physical system to operate in a prescribed 
manner." The set of procedures noted above is as broad as 
the definition. One may choose to perform the design effort in 
the time domain by applying modern state space concepts or by 
direct solution of the system differential equations. Alternatively, 
the complex frequency domain may be selected, employing any 
of the universally accepted procedures such as Nichols and 
Bode plots, root locus, Mitrovic's Method, etc. 

Notwithstanding the achievements which have resulted 
from direct application of the above procedures, it must be 
admitted that they have inherent limitations . With the 
exception of simple systems , time domain analysis virtually 
demands the assistance of a digital computer. In a sense, 
this is a disadvantage because computers are neithe r available 
to all people at all times nor is it economically practical to 
employ a computer for all computation. The complex frequency 
schemes are considered imperfect because they are time 
consuming, i.e. , essentially graphical in nature and often 
based on trial and error, and normally limited to one or two 
variables . 
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In view of the above discussion, it appears that there is 
need for still another design tool. However, in order for a 
new procedure to even merit consideration, it must be 
relatively free from the faults and limitations cited above. 
Specifically, the method should be analytically simple, i.e. 
extend the range of problems which can be solved without 
resorting to graphical techniques or the exclusive use of a 
computer. It should also be rapid, capable of handling several 
variables, specific and accurate. The remainder of this paper 
is devoted to the development and discussion of a design 
procedure which is believed to have these attributes. 
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2 . The Algebraic Method; Concept and Development 

Concept 

The proposed design procedure is based on the premise that 
system design specifications can be transformed into algebraic 
functions of the system parameters . If the specifications are 
independent, N specifications will produce N equations which 
simultaneously constrain the system. Control system design 
would thus revert to the relatively simple problem of solving a 
set of simultaneous algebraic equations. 

The problem of transforming specifications into equations 
is not as broad and difficult as one might initially assume. 

The most common and practical complex frequency specifications 
are; undamped natural frequency (w^^) , darning ratio (3^ ) , steady 
state error (Egg), bandwidth (wj^) , reasonance peak (ivip^) and 
reasonant frequency (wj.) . In many cases, both ^ and Wj^ are 
specified . This is equivalent to specifying the location of a 
complex pole pair and generally implies that these poles should 
dominate the system response. If attention is restricted to this 
situation, then the specifications of primary interest are , 

Wn, W]^ and Egg because Mp^ and Wj. can be at least 
approximately expressed in terms of Jr and w^. The remainder 
of this section is therefore devoted to the development of methods 
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which can be used to obtain system constraint equations from 
specifications for jr , Wj^, Wj^ and Egg. 

Root Location and 

The coefficient equations associated with the specification 
of a complex pole pair can be developed in two ways. The first 
method is based on the root -coefficient relations of an 
algebraic equation; the second is based on Mitrovic's Method. 

Assume the systems' characteristic equation (CE) , obtained 
by setting the c^nominator of the closed loop transfer function 
equal to zero, is the polynomial. 







Sk - O 



where s is a complex variable and defined as 



( 2 - 1 ) 



5 = _ 



( 2 - 2 ) 



and the a|.*s represent the real coefficients which are functions 
of the system parameters. 

Since equation (2-1) is of order n, there must be n roots 
which satisfy the equation. Therefore, an equivalent expression 

i . 

for equation (24i) is 









(2-3) 



where represents the ith root of the polynomial. By expanding 
equation (2-3) and equating the coefficients of like powers of s in 
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equations (2-1) and (2-3) , n equations can be obtained which 
relate the roots and coefficient of the original polynomial. 

If ^ ^ and Wq define the location of the specified complex 
pole pair, it is permissible to set 

(2-4) 

Z. — ~ / y/— ' 

(2-5) 

Substituting these values of and R2 into the n root-coefficient 
relations obtained previously gives n equations in terms of the 
remaining n-2 roots and the coefficients of equation (2-1) . By 
algebraic manipulation, it is then possible to successively 
eliminate each of the remaining n-2 roots and finally obtain 
two equations involving only the coefficients. Since the 
coefficients are explicit functions of the system parameters, 
any combination of the parameters which satisfy these equations 
will generate the desired pole pair. 

Example 2 . 1 

Let the system characteristic equation be 
-h + dtS Qo 

Find the relations which must exist between the coefficients in 
order to locate a pole pair withj^ = 1/2 and Wq = 10. 

Denote the three roots of the CE as R^, R2 and R3. Then 
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and ct , — y- (Rt-*- Rt. ) 

— RtRs./R3 

= R, 'f' Rz. ~*~Rs 



Setting 


- ■tz 


^ J fo / ^/y’ 






- sT 


- J /o / W 




Then 


/oo (Qz-/o) 




(1) 




/O(ar-/0) 


— Oo 


(2) 




/ooRs 


= ao 


(3) 



Equations (1) and (2) are the desired coefficient relations; 
equation (3) provides information about the third pole. 

While the root-coefficient method works well with low 
ordered systems (and even provides information about the 
remaining n-2 roots), it requires an inordinate amount of 
computation when applied to systems of higher order because 
the coefficient relations must be extracted from a set of non- 
linear equations. This difficulty can be overcome (at the expense 
of losing information about the remaining n-2 roots) by using 
Mitrovic's Method. With this method, a kth ordered system can 
be handled as easily as one of third order. The method is fully 
developed in (3) and (5). Only a brief review is given here. 

The argument begins with equations (2-1) and (2-2) defined 
previously. Two independent equations in aj^, w^ and '-jArcbn be 
obtained by substituting equation (2-2) into equation (2-1) and 
setting the real and imaginary portions of the resulting 
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expressions equal to zero. Through additional manipulation and 
the introduction of a new variable, the equations take the form 



4< (?") ~ o 

K-O 


(2-6) 


OZ K (f)~^ 

K"0 


(2-7) 


Where 4^ (f) = - F 4^., (?> (T) ] 


(2-8) 


with 4o(F)^0, 4,(f)=-l 


“ 4< (?) 



Other values of ) for various values of y are tabulated 

in Appendix I. 

Equations (2-6) and (2-7) are known as the generalized 
Mitrovic equations and can be used to obtain the coefficient 
relations required to locate a pole pair as specified by JT q 
and Wq. 

Example 2 . 2 

Work example 2.1 using Mitrovic' s Method. 

For a third order system, the generalized Mitrovic equations 
are 

oj^4^ cT) 4>, (?) ^ ^3 ^ 

Substituting Wj^ = Wq = 10 and using values for jz(j^(l/2) from 
Appendix I, the above equations simplify to 

Qt = / O . * — /OO fOOO 
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The results are equivalent to those obtained previously. 

Although only a single pair of complex poles is being 
considered here, it should be noted that either method is 
capable of handling any number of real or complex poles 
(providing, of course, that the total number considered is less 
or equal than the order of the system). For example, assume 
two complex poles pairs and two real roots of a sixth order 
system are specified. The root -coefficient method will provide 
six constraining equations by direct substitution of the given 
values. If Mitrovic's Method is used, four constraint equations 
dependent oh the complex poles can be obtained by repeated 
use of equations (2-6) and (2-7). The remaining two equations 
which are dependent on the real roots are produced by repeated 
use of 

**y> 

(2-8) 

obtained by letting s = -Q“ in equation (2rl) . 

Steady State Error (Egg) 

There are three different steady state error constants which 
may be specified. They are 

position error constant = = lim G(s) (2-9) 

s-*o 

velocity error constant = = lim sG(s) (2-10) 

s-»o 

acceleration error constant = k^ = lim s^G(s) (2-11) 

s-^o 
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The above definitions are only valid for systems which have all 
their roots in the left-half of the s -plane (because the 
definitions are merely restatement of the final value theorem 
for complex variables) and unity feedback around the forward 
loop transfer function, G(s). To satisfy the last requirement, 
it may be necessary to manipulate the block diagram into an 
equivalent G(s) block with the required unity feedback. 

The constraint equation is obtained by equating the specified 
value of Egg with the algebraic expression resulting from 
application of equation (2-9), (2-10) or (2-11). The procedure 
is demonstrated in the following examples. 

Example 2 . 3 

For the system shown in figure 2-1, obtain an expression 

between K, Z and P for k^ = 6. 

a 

For the given system 







Applying equation (2-11) 






Ki: 
30 R 



which, for kg = 6, reduces to: 18 OP = KZ 

Example 2.4 

For the system shown in figure 2-2, derive an equation in 
terms of Kj^, K 2 and which will give = 8. 
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Since the system does not have unity feedback, insert 
fictitious positive and negative feedback paths around the 
system. Geq(s) can then be found by including the positive 
feedback .branch . 



Applying equation (2-10) 





-t- K, /<J ) •#■ s (s-h K.Kx) - K, 






s <^£0^5) = O 

s-*-o ^ 



Therefore the specification cannot be satisfied. 

Bandwidth (wj^) 

Bandwidth will be defined as the angular frequency, Wj^, 
at which the magnitude of the closed loop transfer function has 
a value of 0.707, i.e. , - 3 db bandwidth. The transformation 
of the specification into the desired constraint equation can be 
accomplished by following the usual procedures: 

A. Find the closed loop transfer function and set s = jw. 

B. Find the magnitude of the resulting expression and 
equate it to 0 . 707. 

C. Substitute w = wj^ and simplify. 

Example 2 . 5 

For the system shown in figure 2-3, derive a relation in 
terms of K and K^. which assures the system has a bandwidth 
of 10 radians/second. 
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Following the prescribed procedure: 



C^s) _ K 

^(s) -h £"s -i-S K 



C:<jta) /cf 



For w = Wj^ = 10 , 



• • 







/O**' —ZOOK 'i- /00(S~hKKr)^ 



t 



The constraint equation obtained in the preceeding example 
for a second order system is typical of the expressions which 
are developed from bandwidth specifications . Since the 
expressions increase in complexity for higher ordered systems 
or for systems having closed loop zeroes, most design problems 
involving bandwidth specifications will usually require the 
assistance of a computer or a graphical technique. (A.n excellent 
graphical procedure based on Mitrovic's Method is given in (ft),) 

If such aids are not available, a convenient, and relatively fast, 
alternative solution is to initially ignore the bandwidth 
specification and design the system to meet the remaining 
specifications. The results can then be inserted into the system 
bandwidth equation to determine if they also satisfy the bandwidth 
requirement. If the specification is not satisfied, the magnitude 
of the deficiency may suggest a more reasonable structure for the 
next attempt. 
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(5+6)(s>f5) 



(s2SL 

(s+p> 



FIGURE 2.1 

BLOCK DIAGRAM FOR EXAMPLE 2.3 




BLOCK DIAGRAM FOR EXAMPLE 2 . 4 



QC5) 




BLOCK DIAGRAM FOR EXAMPLE 2.5 
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3. Application of the Algebraic Method. 

As stated previously, the basis of the proposed design 
method is the transformation of N specifications into a set 
of N simultaneous constraint equations that are explicit functions 
of the system parameters . Since the specifications are 
independent, a specific design will generally require that the 
system have N parameters. The insertion of sufficient 
parameters, usually compensating elements, is an integral part 
of the design procedure. 

It should be noted that the mere selection of sufficient 
parameters does not guarantee a satisfactory solution. A 
converse of this statement is also true. An unsatisfactory 
solution does not imply that a solution does not exist; one may 
possibly be obtained by employing a different set of parameters. 
These statements suggest that there is some trial and error 
associated with the procedure. This is correct. Hbwever, the 
amount of trial and error required can be minimized through 
experience and application of control theory fundamentals. In 
this respect, root locus sketches are particularly valuable in 
obtaining clues to the nature of the required parameters. 

The algebraic method has the inherent capability to handle 
any number of specifications or parameters. Therefore a 
parametric solution can generally be obtained by simply inserting 
more parameters into the system than there are design 
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specifications. If M parameters are selected for use in the 
design of a system subject to N specifications, for M>N, the 
results will be parametric and expressable in terms of any M - N 
parameters . 

Example 3 , 1 

Given G(s) = K/s(s + 4), design* a system which will 
produce a pair of complex poles defined by JT q = 1/2 and w^ = 4 
with = 4 . 

It will be stated without proof that K = 16 satisfies all 
specifications. 

This seemingly trivial example has major significance. 
Specifically, one parameter successfully satisfies three 
specifications which is a direct contradiction of the previous 
discussion. It would appear that adherence to the general rule 
of N parameters for N specifications could possibly result in 
the incorporation of redundant material into the system. This 
does not seem likely because of the following observation: 

In all cases studied, parameters which did not assist in 

the solution tended to be eliminated . 

Examples 3.2 and 3.3 will be used as illustrations of this 
observation. 

Example 3 . 2 

Insert cascade compensation, G^(s) = (s+Z)/(s+P), into 
the system given in example 3.1 and redesign the system to 

*Note: All designs will employ unity feedback between output and 
input (in addition to any other specified feedback) . 
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loot 






/eg y- >C5 

^3 ^ 31 4- ^ PS 



Applying equation (2-10): 



K> 



= -/^. g 



( 1 ) 






System CE: 

5 2 -h^(>+-tP) -4- 5 (f^P-i-K ) = O 

Mitrovic equations: ^ < (2) 

/('S ^ /C>i^-hP)~C^ (3) 

A simultaneous solution of equations (1), (2) and (3) gives: 

K = 16 P = Z 

Since the p61e and zero of the filter are coincident, the filter 
performs no function and can therefore be eliminated. 

Example 3.3 

For G(s) = 50/s(s + 5), it is possible to locate poles at 
-Sij 5 ( j?*Q = 0.707 and = 5 (2) by using either 
tachometer feedback (K^ = 0. 1) or a cascade compensator (Z = 5, 

P = 10) . In each case, = 5. Design a system using both 
tachometer feedback and a series compensator to satisfy the 
same specifications. 

For combined tachometer feedback and cascade compensation: 






•TO 



t-"5^ ^ P •i-S’i-SO Kr) fS (SPf SO Kri) 



Applying equation (2-10); <v = Soz/(SP ^ (D 

21 



System CE: 

-h d^(P’tStSoKT)i-sisP-hSOKT=t tSo)^SO^ 

Mitrovic equations: So 'B - SO SOKt)-‘SOO (2) 

SP fSrOkCrlt -hSO « /O ( R-t S SO Kr)- SO (3) 

Simplifying: Z = 5 (4) 

P = 10(1 - Siq.) (5) 

The result obtained in this example is doubly significant. 

First of all, it indicates that the procedure is capable of 
detecting conditional as well as complete redundancy. In 
example 3.2, the filter was redundant for all values of the 
compensator pole and zero. In this example, however, 
equation (5) shows that the filter is redundant only when 
Kj = 0.1. Similarly, the tachometer may be discarded when P = 10. 
The second point of interest is the fact that a unique solution 
may not exist. In this case, any positive values of P and Kj 
which satisfy equation (5) also satisfy the specifications. Which 
solution is best? At this point, it is a matter of conjecture. 

However, this example will be rediscussed in Section 4 after a 
means for determining the best design has been developed. 

Example 3 . 4 

Given the forward loop transfer function, G(s) =— ; — 77 — — r } 

design series compensation to give roots with f q - 1/2 and 
Wq = 6. The minimum acceptable value of is 4.15. 
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For cascade compensation: 



so {s) = 



K 



S‘^+ 7-SPs^-f (7.SP *■ S.S^Ps 



Applying equation (2-10): Km — ^./S « 



- 



3.SP 



( 1 ) 



System CE: 

S'^ V' 7.S Ps"^ -h ( 7.SR^^.S)-^S(3^ P7-K) * K^-O 



Mitrovic equations: 



Acz- = s^p- 



( 2 ) 



K = Hh^P^ /9S- (3) 

Solving equations (1), (2) and (3) gives: K= 1475 

P = 37.9 
Z = 0.373 

Example ?2 > . 5 

The bandwidth for the system in example 3.4 is about 6.5 
radians/second. Redesign the system so that the bandwidth is 
reduced to 4 radians/second. 

The introduction of another specification requires the 
insertion of another parameter. Use tachometer feedback around 
the original plant and series compensator. Then 






(7. £--rP)-t-s^{^,S-^7S-p-y^y^7)fs(s^Pf-K<T^) 



Applying equation (2-10): Km 






( 1 ) 



System CE: 

(7‘S-fP) -h (^'■S-»-?-SPt-KKT)i-s(3-SPfHf(-T^TK)-hKi »0 
Mitrovic equations: ^ ^2(3.s^y.s-Pi-k*fr)-cHP*7.£) (2) 

3.S Pf-KiCrir = (o{3.^-r y.SPtKKr)- ^ 2 ) 



System bandwidth equation for Wj^ = 4: 



1/2 = 






/6xrs 



[Ki:-iUi3.£fl.SPtlilCr')-ZS(^^->- K~ iG (P ->-7.S)y 



C^) 



Solving equations (1) through (4): 



K = 252 


P = 5.93 


Kt = 0.16 


Z = 1.015 
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4. Sensitivity as a Specification. 



Attention is now directed to the question posed during the 
discussion of example 3.3. If more than one design satisfies 
the specifications, which is the best design? The fact that a 
specific solution was not obtained indicates the need for another 
equation which, for the system being developed, implies the need 
for another specification. Since the design specifications have 
already been incorporated into the solution, it follows that the 
additional specification should be general in nature and 
applicable to all systems. 

One choice for a general control system specification is 
Sensitivity which will be defined as 



to a parameter Q . " (This definition is preferred by the author 
because it imparts more physical meaning than most other 
definitions of sensitivity found in the literature, i.e. , it states 
that sensitivity is essentially the ratio of the percentage change 
in the function to the percentage change in the parameter.) The 
reason for this choice is the fact that a sensitivity specification 
will force the final design to be realistic. Specifically, it requires 
the consideration of changes in the parameters of the physical 
system which will result because of tolerances, aging and 




(4-1) 



P 

where Sq is read as “the sensitivity of a function P with respect 



environment. On this basis, the best design is the one which not 



only satisfies the design specifications, but is also able to 
maintain the desired performance in spite of undesired parameter 
variations , 

Root Sensitivity 

Use of equation (4-1) requires definition of the function P. 
Since sensitivity is being introduced as a means of insuring 
compliance with the design specifications, P should ideally 
be a function of all the design specifications. However, in 
this presentation, only the root specification, i.e. , jS’ q and 
Wq, will be considered. Hence, the discussion will refer to 
the root sensitivity of the system. (Hereafter, root sensitivity 
and pole sensitivity will be used interchangeably.) 

Restricting attention to root sensitivity does not completely 
preclude the determination or discussion of a best design. In 
fact, there are two arguments which partially justify (and might 
even suggest) this simplification. One argument is that root 
specifications are more stringent than either the bandwidth or 
error coefficient specifications. The latter two are frequently 
specified in terms of inequalities, i.e. , greater or smaller than 
a certain minimum or maximum value. Therefore, some tolerance 
can be built into the design by working on the "safe side" of the 
specifications. Another argument is based on the root- 
coefficient relations of the system characteristic equation. If 
the mobility of two of the systems' closed loop poles are 
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restricted by means of a sensitivity constraint, then the mobility 
of the remaining closed loop poles is also restricted. This 
statement stems from the fact that the specified poles have 
corresponding roots which must be factors of the characteristic 
equation. Now, in order for these roots to be perfect factors of 
the characteristic equation, the coefficients of the characteristic 
equation must have remained reasonably constant for small 
parameter changes. This implies that the remaining roots, and 
hence the remaining poles, also remained reasonably constant. 
Since both bandwidth and the error coefficients are partially 
dependent (and, in fact, totally dependent for low order systems 
without zeroes) on the characteristic equation, it follows that 
their sensitivity functions bear an approximate relation to the 
root sensitivity functions . 

Root Sensitivity Equations 

Although the practical importance of sensitivity has long 
been recognized, it has only been in recent years that effort 
was expended in the development of sensitivity as a useful tool. 

I I 

The work of Kokotovic and Siljak, (2) , is considered to be a 
major advancement towards this goal. In essence, they used 
the definition of sensitivity, equation (4-1), in conjunction with 
the generalized Mitrovic equations, equations (2-6) and (2-7), 
to obtain algebraic expressions for (what the author considers 
to be) the components of root sensitivity. Their results are not 
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only particularly suited for use in the design of control systems 
based on a specified pole pair, but also simplify sensitivity 
computations because they permit the exclusive use of real 
variables in calculating the sensitivity of a function containing 
complex variables. 



Since the Kokotovic and Siljak sensitivity equations will 
be used in subsequent examples, a brief review of their develop- 
ment will now be presented. Differentiating the generalized 
Mitrovic equations, equations (2-6) and (2-7), with respect to 
a parameter Q gives 



j^k-i( P* ) respectively and are related by the general equation 



Numerical values for for various values of jp ^re given in 
Appendix II. Manipulation of equations (4-2) and (4-3) and 
application of equation (4-1) gives 




(4-2) 




4k(f) = -z[t 



K-> 

<- 2 . 



4>k-x. (T)J 



(4-4) 



Sq ^ ^ Si — Q Cl 



(4-5) 




(4-6) 
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where 






' o 




(4-7) 


B, = 4*K-i Cf) 

o 


o 


(4-8) 






(4-9) 



Solving equations (4-5) and (4-6) for Sq and S^n gives 



where 





- Q 




(4-10) 






“5“ 




O ( 3 ^ 




D 


(4-111 



0 = 


A, 


3, 








Bz 


(4-12) 




C. 

Cz 


&i 1 


(4-13) 


Dp = 


A, 




(4-14) 








1 



Equations (4-10) and (4-11) are the desired sensitivity design 
equations. (Hereafter, Sq and Sgn will be referred to as the 
"component sensitivity functions. ") At this time, their true 
simplicity is masked by their dependence on the solution of a 
number of auxilliary equations. Suffice to state that the solutions 
of the auxilliary equations become exceedingly simple when the 
location of a pole pair is specified ( ^ q and Wq) and the 
coefficients of the characteristic equation are known (a^) . This 
will be clearly demonstrated in subsequent examples. 
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Intrepretation of Sensitivity Values 
Equations (4-10) and (4-11) may be used in two ways. 
Substitution of a given set of parameter values into the right 
hand side of the equations gives a definite sensitivity value 
for a specific set of parameter values. Alternatively, the 
assignment of a sensitivity value to the component sensitivity 
functions give the parameter relations which must be satisfied in 
order to obtain a system of the specified sensitivity. 

In either case, the final design will be dependent upon the 
accepted/assigned value to sensitivity. It follows that some 
guidelines for reasonable sensitivity values must be established. 
Since the literature is lacking in this area, the following guides 
(which are consistent with the original interpretation of equation 
(4-1)) are proposed by the author: 

A. Signs (+ or -) of sensitivity values indicate sense 
only. Therefore, the minimum value of sensitivity is zero, 
implying that a function is totally insensitive to variations of a 
parameter. 



B. js^ = 1 indicates aUinear variation in P due to changes 
in Q.i.e. , a 5% change in Q results in a 5% change in P. This 
is considered to be the maximum permissible value of sensitivity; 
if possible, systems should have sensitivity values much less 
than one. 
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c. 



S ^ >1 indicates amplification of the parameter 

Q I 

variation. This is generally considered undesireable for a 
control system (but may have application in some detection 
schemes) . 

D. As indicated in A above, sensitivity values have vector 
properties. Therefore, the total root sensitivity, is 

appropriately expressed as 

(S^T)2 = (s5>2 + (S}^n)2 (4-15) 

Sensitivity as a Specification 

The concept of sensitivity has been developed to the extent 
that it can now be used as an additional system constraint 
subject to analysis by the algebraic design method. The 
remainder of this section will be devoted to the design and 
discussion of specific systems which must satisfy sensitivity, 
as well as other, specifications. 

Example 4 . 1 

Reconsider example 3.3 in view of a specification requiring 
minimum sensitivity of the specified poles with respect to 
variations in system gain, K. 

Considering the system gain, K, as a parameter, the system 
characteristic equation is 

( 1 ) 

The sensitivity equations will now be devloped in detail to 
demonstrate the solution of the auxilliary equations. Applying 
equations (4-7) and (4-8) gives 
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( 2 ) 



/?/ = / Soo — /OO 

/?2 ^ Q, — 'pso\fl 

Bt = ^Got/2 03 

Sx. - /OO < 72—2 OCO 

after substituting = Wq = 5 (2)^/2 and values for j^j^(0.707) 
and 707) from Appendices I and II. Although specific 

values of P and Kj. are unknown, inspection of equation (1) in 
conjunction with the results of example 3.3 show that the 
coefficients of the characteristic equation must have the 
following values for all solutions 

<5 ® ■=■ Z£0 <!, = /OQ • !/ ~ ! 

•« t ‘* i V 

which, when substituted into equation (2), gives 



/?/ 



At.^zso\fl 



.-£00 



(3) 



Taking the partial of each coefficient in equation (1) with 



respect to K gives 

<s><7< 



- o 






(4) 



Applying equations (4-9) and substituting, as above, for w^ 
and (0.707) gives 

C, = 2 - 50 kr <=z~ So/TKt 1 1 ) 

which, for Z = 5 as obtained in example 3.3, reduces to 

C, ■=. S‘Ci-fOKr) Cz - .C(/r -/) (5) 

Applying equation (4-12) and substituting equation (3) gives 

D ■= -ZSO X /ci^ (6) 

Applying equation (4-13) and substituting from equations (3) 
and (5) gives 
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( 7 ) 



Applying equation (4-14) and substituting from equations (3) 

and (5) gives x 

:=^-iZ^c/T (f-foK.r) ( 3 ) 

Substituting Q = K= 50,/’ = JTq = 0.707 and equations (6) , 

(7) and (8) into equations (4-10) and (4-11) gives 

sS^ = (/oKr — f ) 

‘/z (10) 

Equations (9) and (10) are the desired design equations. 
Since only one additional constraint equation was needed, they 
would normally be squared and added to obtain an expression 
for the total root sensitivity. In this case, however, S^nis 
observed to be a constant. Therefore, the minimum pole 
sensitivity with respect to K can be found by setting equation 
(10) equal to zero. This gives = 0.1 which corresponds to 
the solution of tachometer feedback only. 

For K|. = 0.1, S]^n = 1/2 and Sj, = 0. Therefore, the analysis 
predicts that a ten per cent change in K should produce a five 



per cent change in w^^ while ^ remains constant. These 
predictions can be mathematically verified. The system 
characteristic equation for the specified value of is 





+ s (K Kt -i- S) -r K. 


= o 


For K = 50: 


Wq = 7.07 


f 0=0-707 


For K = 55: 


w = 7,42 
n 


Jf = 0.708 
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Parameter Selection 



The preceeding example raises several questions which 
have yet to be discussed. The first deals with the selection of 
a parameter on which to base the sensitivity analysis. One 
choice might be the parameter which is considered to hav^ the 
widest variations during expected conditions of operation, e.g. , 
the gain of an electronic amplifier. Another choice might be 
similar to the first but restricted to plant parameters. This 
assumes the designer is able to select and/or control the 
sensitivity of the required compensating elements. There is no 
guarantee that either is correct. It may be possible that the 
system is naturally insensitive to the parameter likely to have 

I 

the widest variations in value while it could be highly sensitive 
to a parameter that is expected to have only minute variations 
in value. It is also possible that, by minimizing the sensitivity 
due to one parameter, the effect of another isiatuplified. . jThe 
answer, then, requires investigation of all possible parameters. 
It is believed that a thorough investigation of a particular type 
of system could lead to the establishment of general guides for 
the selection of the appropriate parameter(s) in other similar 
systems. Since such an investigation is not consistent with the 
purpose of this paper, the author will be content to arbitrarily 
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select any parameter which will facilitate demonstration of the 
algebraic method. 



It should also be noted that there are no restrictions as to 
the number of parameters which may be similtaneously considered. 
If the sensitivity of pole P with respect to parameters Qj , Q 2 
and Qg is desired, the component sensitivity functions for 
each parameter may be individually calculated . The total 
component sensitivity function due to all three parameters is 
the algebraic sum of the individual components, e.g. , 



QxQs 









<^3 



. Since the sense 

of each component must be considered, there is a possibility 
that individual components may cancel. Another fortunate 
circumstance is that the computations for a situation considering 
the root sensitivity with respect to N parameters do not increase 
N-fold. This is because equations (4-7), (4-8) and (4-12), i.e., 
Aj, A 2 , B^, B 2 and D, remain the same for all parameters. 

Minimum Sensitivity 

Another point of interest in example 4.1 is the ease with 
which the condition of minimum sensitivity was obtained. This 
is not a unique situation because it was possible to obtain 
results by inspection in other systems studied also. However, 
it seems quite likely that for many systems the parameter 
relations for minimum sensitivity must be obtained by setting 
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the derivative of the sensitivity function equal to zero. (Note 
the problem if the sensitivity with respect to more than one 
parameter is being considered.) In such cases, care must be 
taken to insure that the function to be differentiated is in a 
form which will indicate minima near zero. This is necessary 
because minimum sensitivity was defined to be zero. The above 
stipulation can be conveniently accomplished by considering 
either the absolute value or the square of the component 
sensitivity function, i.e., or vice ^ 

Use of the total sensitivity function, equation (4-15), is also 
appropriate . 

Discussions thus far have repeatedly used the term "minimum 
sensitivity." This does not necessarily mean that minimum 
sensitivity must be attained by, or even an objective of, every 
design. The specified sensitivity value, whether it be zero or 
two, should be consistent with the other specifications. For 
example, if the pole locations are only required to be accurate 
within five per cent for satisfactory system response and the 
parameter of interest is also expected to have five per cent 
variations, a sensitivity value of one is sufficient. Of course, 
if analysis shows that a lower value is attainable, the design 
should be based on the lower value provided it can be achieved 
without additional cost or degradation of some other aspect of 
system performance. 
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The nature of the sensitivity function in the vicinity of the 
design value is also of concern. If the design value is 
associated with a point on the sensitivity function curve where 
the curve is very steep or irregular, the system is likely to be 
more sensitive than an alternative design based on a (reasonably) 
higher value which corresponds to a smooth portion of the 
sensitivity function curve. 

Example 4.2 

A portion of a control system has the forward transfer 
function, G(s) = K/s(s + 2) (s +4). Dominant closed loop poles 
with J- o “ 1/2 ^nd Wq = 5 are required for the desired response. 
Overshoot is of prime importance; therefore , ^ ^ should be 
maintained within 1% accuracy for ’a- anticipated gain variations 
of 5%. Velocity and acceleration signals are required elsewhere 
in the system and are therefore available for use as feedback 
signals. 

The problem as stated gives 3 specifications (pole pair 
and S)t) and 3 parameters . (Note that an additional specification 
concerning S]^n would normally require the insertion of an 
additional parameter whereas a specification regarding the total 
pole sensitivity would not. Also note that the specifications 
permit Wq to vary, if necessary, in order to maintain ^ 
approximately constant.) 
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System CE: -S* (i> K Ka.) -t s KKt) +• K = o 



Mitrovic equations: 



2S t- ZS KKa 



KKt ^ zz r s 



( 1 ) 

( 2 ) 



Applying equations (4-7) thru (4-9): 



— 37S— Soox 
&, = ZS-O 
C"/ - !~ZS /< 



/? a = SOo2, — ^o, 
Sa. -Sd)aa.-^00 

C a. * ~ ) 



(3) 



Unlike example 4.1, the coefficients of the characteristic 
equation do not have fixed numerical values. Therefore, the 
solution of the determinants for D and Dj^, equations (4-12) and 
(4-14), will be hampered by the presence of all 3 parameters. 

At this point it will be convenient to solve equations (1) and (2) 
for K and K^. Then equation (3) can be expressed in terms of K 
only. 



Kn =• 

Therefore 

/), = tZS-ZK 



\<-z.cr 

ZS‘K 



I/- = 

5K 



(4) 



^ K I^S 



B, ^ 2 SO ^ZK-ZSO 

C. = 2S/k Cz= —U0/< 

T 

Using equation (5) to solve for D and D|^: 

O^-ii-K^+SooK-ZSO^ ^(i9SK-I5SxIZ5) 

Using equation (6) to solve equation (4-13: 

- 97 . 5 - (K-BC.S) 

K^— IZStC -r IZS^ 



(5) 



( 6 ) 



(7) 



Inspection of equation (7) shows that s|^ = 0 for K = 86.5. Since 
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the only requirement is that Sj^ ^ 0.2, this gain setting should 
be considered as a possible design value. For K = 86.5, 

= 0.397 and = 0.0284. For these values, the third closed 
loop pole is located at 'CT“= -3.4. Therefore, the complex poles 
are not completely dominant. In addition, figure 4-1 shows the 
sensitivity function is very steep at the point of minimum 
sensitivity. Again using figure 4-1, the system will be designed 
on the basis of K = 600. The other design values are; K^- = 0. 238 
and = 0.0384. The complex poles are assured of being 
dominant because the third pole is now located at CT"= -24. It 
is stated without proof that S^n also approaches zero as K 
becomes very large. It is also less sensitive than sj^, having 
a value < 0.05 at the design gain setting. 
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l5^l 




s £/js/7~/ \// ry (S%) vs (k) 

/-OR B •4—Z. 
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5 . Conclusions 



The purpose of this paper was to develop a design technique 
based on the premise that design can be accomplished by solving 
a set of algebraic simultaneous equations derived from the 
system design specifications. In addition to demonstrating the 
feasibility of this concept, the examples in the preceeding 
sections also showed that the method developed has certain 
desirable characteristics. Specifically, the algebraic method 
is analytically simple, capable of handling any number of 
parameters, rapid, specific and accurate. 

The method is in no way restricted to the exclusive use of 
the design specifications considered in this paper. The 
specifications considered were simply selected to enhance the 
presentation. In general, a system can be designed to satisfy 
any set of independent specifications providing the specifications 
can be expressed as algebraic equations. It is obvious, of 
course, that the solution of the resulting equations becomes 
more difficult as the number of specifications increase. For this 
reason, it is suggested that analysis be normally limited to 
specifications which are considered to be prime, i.e. , those 
which are exactly specified as opposed to others stated in terms 
of inequalities or approximations. The results of this analysis 
can then be used to determine if the secondary specifications are 
also satisfied. 
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The inclusion of sensitivity equations into the set of system 
constraint equations was considered especially significant. Here 
again, the use of sensitivity equations is in no way dependent 
upon the assumptions made in Section 4 concerning the relation 
of root sensitivity to system sensitivity, i.e. , use of root 
sensitivity as a performance index. This was logically 
demonstrated in example 4.2 where a root sensitivity component 
was of particular concern while the system sensitivity had no 
restrictions whatsoever. Although only root sensitivity was 
discussed, there is no apparent reason why similar restrictions 
cannot be imposed on other design specifications. 

As previously indicated in Section 4 , there is a definite need 
for the study of specification/parameter sensitivity relationships 
before sensitivity can become a useful design tool. At the 
present, the sensitivity effects of all possible parameters must 
be evaluated with respect to each specification and to each other 
before a meaningful design can be achieved. Analysis of even 
a simple system on this basis would require a prohibitive amount 
of computation because of the large possible number of parameter/ 
specification combinations available. 

"Seek Simplicity and Distrust It" 

Alfred North Whitehead 
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